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ABSTRACT 


This  report  contains  a  calculation  of  breakdown  voltage  for 

various  types  of  junctions  using  the  most  recent  and  probably  the  most 

accurate  measurements  available  for  electron  and  hole  ionization  rates. 

In  addition,  there  is  a  calculation  of  multiplication  at  voltages  less 

than  the  breakdown  voltage.  This  calculation  is  useful  to  obtain  the 

BV  for  n-p-n  and  p-n-p  transistors.  It  is  found  that  for 
CE 

similar  collectors,  the  BV  for  an  n-p-n  transistor  is  about  half 

CE 

the  value  for  a  p-n-p  transistor. 

In  addition  to  the  multiplication  calculations,  there  is  included 
a  discussion  of  previous  work  done  on  this  contract,  and  an  attempt  to 
reconcile  the  discrepancy  between  the  results  of  calculation  of  electron 
temperature  at  high  and  low  fields.  Inclusion  of  transverse  acoustical 
phonon  scattering  results  in  a  better  fit  of  experiment  and  theory  than 
just  use  of  optical  phonons  and  longitudinal  acoustical  phonons  but 
there  is  still  a  significant  discrepancy. 
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SYMBOLS  FOR  APPENDIX 


electron  mass 


m 

o 

m^  =  transverse  electron  mass 

N  =  total  number  of  electrons  in  the  conduction  band 
N  =  number  of  transverse  acoustic  phonons  on  the  zone  edge 

ci 

t  h 

=  number  of  phonons  of  i  mode 
Nq  =  number  of  optical  phonon  at  zone  center 

n^  =  number  of  ions  per  unit  volume 

pi(c°s  0)  th 

=  Legendre  polynomial  of  i  order 

— > 

p  =  electron  momentum 

Pq  =  momentum  associated  with  electron  drift 

q  =  phonon  momentum 

qQ  =  greatest  phonon  momentum 

C2 

r  =  -  ,  dimensionless  constant  of  the  order  of  a  few  times  10-2 

18D2 

e2 

R,  =  -  ,  adjustable  parameter  1 

1  4r62 

o 

e2 

R„  =  -  ,  adjustable  parameter  2 

2  4r62 

a 

s  =  longitudinal  sound  velocity  in  the  crystal 

T  =  effective  electron  temperature 

Tq  =  lattice  temperature 

At  =  incremental  change  of  electron  temperature 

V  =  crystal  volume 

v  =  drift  velocity  of  the  electron 
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=  Debye  temperature  of  low  energy  acoustic  phonons  near  zone 
center 

=  Debye  temperature  of  transverse  acoustic  phonon  near  zone  edge 

th 

=  Debye  temperature  of  i  mode  of  phonon 
=  Debye  temperature  of  optical  phonons  at  zone  center 
=  Debye  temperature  of  longitudinal  optical  phonon  near  zone  edge 
=  effective  electron  mobility 

=  electron  mobility  due  to  low  energy  acoustic  phonon 

=  electron  mobility  due  to  transverse  acoustic  phonon  near  zone 
edge 

=  electron  mobility  due  to  i  mode  of  phonon 
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^aA 


cA 
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=  electron  mobility  due  to  optical  phonon  at  zone  center 

=  zero  field  electron  mobility  due  to  low  energy  acoustic  phonon 

=  zero  field  electron  mobility  due  to  transverse  acoustic  phonon 
near  zone  edge 

=  zero  field  effective  electron  mobility 

=  zero  field  electron  mobility  due  to  optical  phonon  at  zone 
center 

fi 
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I,  INTRODUCTION 


This  report  is  a  final  report  on  contract  DA36(039)SC-85339.  Some 
practical  technological  applications  of  studies  of  ionization  rates  in 
silicon  are  described:  l)  prediction  of  breakdown  voltage  and  capaci¬ 
tance,  and  junction  width  at  breakdown  for  a  variety  of  junctions,  and 
2)  calculation  of  electron  and  hole  multiplication  rates  at  voltages 
less  than  the  breakdown  voltage.  These  calculations  are  useful  in 
designing  p-n  junctions  and  junction  transistors.  A  summary  of  other 
research  that  has  been  carried  out  on  this  contract  is  also  included. 
The  results  of  some  of  the  studies  have  been  published — in  some  cases 
they  are  in  previous  status  reports.  The  significance  of  the  work  is 
discussed  but  the  reporting  is  not  repeated  in  detail.  Also  presented 
is  a  fairly  extensive  discussion  (Appendix  A)  of  a  theoretical  diffi¬ 
culty  in  relating  the  low-,  intermediate-,  and  high-field  behavior  of 
charge  carriers  in  silicon  to  each  other.  This  problem  is  not  resolved 
but  we  have  attempted  to  clarify  it. 
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II.  BREAKDOWN  VOLTAGE  OF  JUNCTIONS 


One  of  the  most  practical  objectives  of  measurements  of  ionization 
rates  of  electrons  and  holes  in  silicon  is  the  prediction  of  breakdown 
voltage  in  various  junctions.  By  breakdown  voltage,  we  mean  that 
voltage  at  which  the  junction  will  conduct  a  large  amount  of  current 
with  a  trivial  change  in  voltage.  In  silicon  p-n  junctions,  there 
are  enough  data  to  calculate  this  voltage  in  a  "perfect"  junction.  By 
a  "perfect"  junction  we  mean  one  that  has  no  imperfections  such  as  dis¬ 
locations,  precipitates,  interstitial  atoms,  or  vacancies  in  the  high- 
field  region.  A  perfect  junction  as  defined  does  not  exist,  but  there 
are  surprisingly  large  numbers  of  commercial  p-n  junctions  that 
approximate  the  perfect  junction  in  breakdown  properties.  The  result 
of  intentionally  introduced  imperfections  is  generally  to  lower  the 
breakdown  field  of  a  given  p-n  junction,  resulting  in  some  cases  in  a 
lower  voltage  and  in  some  cases  in  a  higher  voltage.  The  theoretical 
possibility  that  the  introduction  of  imperfections  will  increase  the 
breakdown  field  cannot  be  entirely  discounted. 

Junction  multiplication  is  defined  as  the  ratio  of  total  junction 
current  to  incident  current.  Thus,  if  there  is  an  incident  current  of 
electrons  only, 


M  =  i/l  .  (l) 

n  '  n  ' 

Or,  for  an  incident  hole  current 

M  =  I/I  .  (2) 

Breakdown  occurs  when  the  multiplication  values  are  infinite.  By  this 
definition,  of  course,  a  junction  never  reaches  breakdown,  but  for 
practical  purposes  a  breakdown  voltage  can  be  defined.  The  multiplica¬ 
tion  values  reach  infinity  simultaneously;  they  are  given  by  [Ref.  1] 


where 


M  ,  M 
n  p 


a  ,  a 
n’  p 


electron,  hole  multiplication 
junction  width 

ionization  rate  for  electrons,  holes. 


Junction  breakdown  occurs  at 


w 

r* 

r*  w  n 

J  an  exp 

)  <V  “n*1”1' 

O 

_x  _ 

dx  =  1. 


(5) 


It  is  not  possible  to  calculate  the  voltage  at  which  Eq.  (5)  is  satis¬ 
fied  for  arbitrary  variation  of  Q!  a  without  the  aid  of  a  computer. 

P 1  ^ 

It  is  possible,  however,  to  obtain  a  number  that  is  within  a  few  percent 

of  the  right  voltage  for  many  different  types  of  junctions  by  using 

suitable  approximations.  It  is  found,  for  example,  that  in  silicon, 

(X  «  a  .  If  a  constant  ratio  is  assumed  for  a  /a  ,  we  obtain  the 
p  n  p'  n 

condition 


s 


a 


n 


exp 


i)a  dx1 
'  n 


dx 


1 


(6) 


o 


where  Ctp/0!n  =  7  which  is  integrable  for  constant 
dition  for  breakdown  is 

w 


o 


In  7 

7  ~  1 


7.  The  resulting  con- 


(7) 
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Figure  1  is  a  plot  of  In  7/7- 1  vs  y.  This  function  is  a  measure 
of  the  number  of  electron-hole  pairs  that  an  electron  must  produce  in 
an  average  transit  across  the  junction  to  maintain  the  breakdown  dis¬ 
charge.  For  the  400/l  range  of  y  there  is  only  a  lO/l  range  in  the 
required  number  of  electron-hole  pairs  produced  by  an  average  electron 
so  that  the  actual  value  of  breakdown  voltage  is  fairly  insensitive  to 
the  7-parameter. 


FIG.  1.  ELECTRON  IONIZATIONS 
PER  TRANSIT  AT  BREAKDOWN 

(in  7/7-1)  vs  7  (=  ap/an) 
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Figure  2  shows  measured  values  of  ionization  rate  vs  electric  field 
for  both  electrons  and  holes  [Ref.  2],  It  is  apparent  from  this  figure 
that  y  is  a  rapidly  varying  function  of  electric  field.  The  variation 
of  y  with  electric  field  is  not  easily  included  in  the  breakdown  rela¬ 
tion  [Eq.  (6)]  without  the  aid  of  a  computer.  An  answer  that  is 
perfectly  suitable  for  most  purposes  can  be  obtained  for  a  variety  of 
junctions  by  taking  reasonable  limits  on  y  in  Eq.  (7)  and  obtaining 
an  upper  and  lower  bound  on  the  breakdown  voltage.  If  the  upper  and 
lower  limits  are  within  a  few  percent  of  each  other,  the  answer  can  be 
considered  to  be  satisfactory. 


FIG.  2.  IONIZATION  COEFFICIENTS  a  ,  a  ,  AND 

P  ® 

y  (=  a  /a  )  VS  RECIPROCAL  FIELD.  THE  DASHED 

'  p'  n' 

PART  OF  THE  y  CURVE  IS  EXTRAPOLATED. 


5 


SEL-64-009 


The  electron  ionization  rate  is  reasonably  well  approximated  by 


a  =  a  e 
n  oo 


-b/E 


6  _i 

a  =  1.5x10  cm 

00 


b  ss  1.6x10  v/cm. 


(8) 


It  is  relatively  difficult  to  integrate  the  form  of  (8)  in  even  such 
relatively  simple  junctions  as  graded  and  step  junctions.  These  inte¬ 
grals  have  nevertheless  been  calculated  and  are  [Ref.  3] 

w 

\  a  dx  =  a  W  „ .  ,  (9) 

J  n  max  eff  ' 

o 

where  0C  is  the  ionization  rate  at  the  maximum  field  and  W  „„  has 
max  eff 

been  calculated  as  shown  in  Pig.  3. 


FIG,  3.  PLOT  OF  W  ,,/W  VS  b/E  . 

eff'  '  m 
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The  procedure  for  calculating  limits  on  breakdown  voltage  then  is 
as  follows,  For  example,  consider  the  graded  junction:  choose  b/Em  = 
4.0,  Then,  from  Fig.  3, 

=  0.4  and  ^=2.5xl0-6. 

m 

From  Fig.  2, 

OS  =  2.8xl04  and  y(E  )  =  0.25. 
max  ' '  max 


From  Fig.  1, 


amWeff 


1.8 


W  =  6.4x10  ' 

ell 

W  =  1.65x10 


V  =  2/3  E  W  =  45  v. 

'  m 

The  gradient  for  a  junction  that  has  a  depletion-layer  width  of 
-4 

1.65x10  cm  at  45  v  is  obtained  from  standard  p-n  junction  theory 

20  -4  , 

and  is  8x10  cm  .  The  electric  field  at  a  distance  ±W  „„/2  from 

et  f 1 

the  maximum  field  point  is  about  85  percent  of  the  maximum  field,  at 
which  7  is  about  0.13.  If  this  value  of  y  is  used,  we  require 

2.3 

8.2xl0-5 
-4 

2,05x10  cm 


Veff 


W 


eff 


W 


V  =  2/3  E  W  =  55  v 

'  m 

a  =  4.9xl020. 
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The  result  of  this  calculation  using  a  smaller  value  of  7  is  that  the 
breakdown  voltage  for  a  given  gradient  should  increase  slightly.  The 
use  of  the  smaller  value  of  y  is  actually  a  slight  over-estimate  of 
the  breakdown  voltage  since  most  of  the  ionization  occurs  in  the 
distance  W  ^  near  the  high-field  region.  Since  the  procedure  that 
we  are  using  here  results  in  different  gradient  and  breakdown  voltages, 
it  is  not  possible  to  estimate  the  error  introduced  into  the  calculation 
of  breakdown  voltage  by  the  assumption  of  constant  y  from  a  single 
computation.  The  calculations  of  breakdown  voltage  vs  gradient  or 

impurity  density  are  within  about  1  percent  of  each  other  as  determined 
by  numerous  computations  and  using  the  two  extreme  values  of  y.  We 
must  conclude  then  that  the  constant-7  assumption  does  not  introduce 
significant  error  into  the  calculation  of  junction  breakdown.  Figure  4 
shows  breakdown  voltage  as  calculated  by  the  above  procedure  for  graded 
junctions.  Figure  5  shows  breakdown  for  step  junctions  and  Fig.  6  shows 
the  width  at  breakdown  for  graded,  step,  and  PIN  junctions.  Use  of  the 
upper  and  lower  bounds  for  7  result  in  a  variance  that  is  less  than  2 
percent  for  the  entire  range  of  breakdown  that  is  shown.  Figure  7  shows 
capacitance  per  unit  area  at  breakdown  for  graded,  step,  and  PIN 
junct ions . 
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!0JO  102’  1022  lo23 

•  (cm-4) 


FIG.  4.  BREAKDOWN  VOLTAGE  VS  IMPURITY  GRADIENT  FOR  LINEARLY 
GRADED  JUNCTION. 


FIG.  5.  BREAKDOWN  VOLTAGE  VS  IMPURITY  DENSITY  FOR  STEP 
JUNCTION. 
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FIG.  6.  JUNCTION  WIDTH  AT  BREAKDOWN  FOR  A 
VARIETY  OF  IMPURITY  PROFILES. 


FIG.  7.  CAPACITANCE  PER  UNIT  AREA  AT  BREAKDOWN 
FOR  A  VARIETY  OF  IMPURITY  PROFILES. 
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III.  MULTIPLICATION  IN  TRANSISTOR  COLLECTORS 


At  least  as  important  as  the  breakdown  voltage  of  a  p-n  junction 
is  the  amount  of  multiplication  at  voltages  less  than  the  breakdown 
voltage . 

The  switching  response  of  a  transistor  is  adversely  affected  if  the 
collector  voltage  is  sufficient  to  cause  multiplication,  and  in  some 
cases  the  transistor  may  become  unstable  if  the  product  OMM  >  1  where 
aQ  is  the  low-voltage,  small-signal  alpha  and  M  is  the  collector- 
multiplication  factor.  The  maximum  collector-emitter  breakdown  voltage 
is  determined  by  the  voltage  at  which  aQM  =  1 

The  multiplication  at  voltages  less  than  the  breakdown  voltage  can 
be  calculated  as  a  simple  extension  of  the  breakdown  voltage.  In  the 
example  of  the  breakdown  calculation  where  the  voltage  at  breakdown  was 
45  v,  the  width  was  1.65  microns.  The  width  at  another  voltage  can  be 
obtained  from  the  fact  that  the  width  of  a  graded  junction  varies  as  the 
l/3  power  of  voltage.  Thus,  at  30  v,  for  example, 


0 

CO 

— 

1.44  n 

l/E 
'  max 

= 

3.2x10 

w  ... 

/w 

0.35 

eff  ‘ 

4 

a 

= 

10 

max 

a 

w  ^ 

0.5 

max 

eff 

7 

= 

0.08 

1  -  r  s  i1  "  eXP  ["  (1  '  7)amax  weff]) 

n  y 

S  0.36. 
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It  y  ,  is  used  where  7  .  is  the  ratio  of  cc  lo.  at  a 
/min  'min  p'  n 

distance  W  2  from  the  center  of  the  junction,  we  obtain 
eir 

1-  (i/Mr)  =  0,40  so  that  there  is  a  10-percent  uncertainty  in  this 
calculation.  To  obtain  the  hole  multiplication,  it  is  most  accurate  to 
replace  the  electron-ionization-rate  integral  by  one  for  holes.  Thus 


1 


JL  s  1 

Mp  -  (l/7)  “  1 


W 

P 


where  W  ..  is  the  effective  width  for  hole  ionization.  As  a  rough 
p  eft 

approximation,  we  replace  the  hole  ionization  rate  by 

ap  =  108  exp  [ -3. 7x106/e]  . 

This  is  a  reasonable  approximation  for  reciprocal  field  greater  than 
-6 

2.7x10  and  includes  all  junctions  that  have  breakdown  in  excess  of 

about  30  v.  We  do  not  have  measurements  of  hole-ionization  rate  for 

electric  field  less  than  2.7x10  v/cm  (l/E  >  3.7x10  )  at  which  field 

the  ionization  rate  is  about  100/cm.  The  ratio  7  was  extrapolated 

linearly  to  the  low  fields  (i.e.,  less  than  2.7xl05  v/cm)  as  a  basis  for 

the  calculations  of  electron-multiplication  rate.  The  calculation  of 

hole  multiplication  is  similar  to  the  electron  calculation  except  that 

there  is  somewhat  more  uncertainty  in  the  result.  One  expects  that  at 

very  low  values  of  multiplication,  the  value  of  1  -  l/M  should 

n ,  p 

converge  to 

w 

\  a  dx  . 

J  n,p 

o 

The  large  discrepancy  between  a and  a  requires  that  1  -  l/Mn  be 

less  than  about  0.2  and  1  -  l/M  be  less  than  0.01  for  reasonable 

P 

equality  to  exist  between  the  ionization  integral  and  the  multiplication 

functions.  Figures  8  and  9  show  1  -  l/M  and  1  -  l/M  as  a  function 

'  n  p 

of  V/Vg  for  several  breakdown  values.  As  expected,  the  electron 
multiplication  remains  significant  to  much  lower  values  of  voltage  than 
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MULTIPLICATION  AS  A  FUNCTION 


does  hole  multiplication.  These  figures  can  be  used  to  calculate  the 
collector-emitter  breakdown  voltage  BV^jj,  for  n-p-n  and  p-n-p 
transistors,  respectively.  The  value  of  collector-emitter  breakdown 
voltage  bvce  is  obtained  by  finding  the  ratio  V/Vg  =  BVCE/BVCg  that 
corresponds  to  l-l/M=l-a  for  the  transistor  in  question.  Thus 
a  75-v  collector-base  breakdown  transistor  with  1  -  a  =  0.02  will  have 
BV„,/BV-,d  =  0.38  if  it  is  n-p-n  and  BV__/BV__  =0.70  if  it  is  p-n-p. 
The  slow  variation  of  the  curves  of  Figs.  8  and  9  with  breakdown  voltage 
allows  fairly  wide  application. 

The  fast  variation  of  hole  multiplication  near  the  breakdown 
voltage  also  means  that  the  uncertainty  in  determining  the  hole  multi¬ 
plication  from  the  ionization  rate  is  not  a  great  disadvantage  in  the 
application  of  these  curves  to  transistors.  The  mathematical  simplicity 
that  is  gained  from  using  the  approximation  of  constant  y  seems  well 
warranted  in  light  of  these  results. 
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IV.  OTHER  RESEARCH 


A.  MEASUREMENTS  OF  IONIZATION  RATES 

We  have  used  both  photo  generation  as  well  as  electric  emission  of 
electrons  and  holes  into  silicon  p-n  junctions  to  measure  ionization 
rates.  The  results  of  the  ionization  measurements  are  given  in  Fig.  2. 

It  is  not  possible,  in  the  arrangements  we  are  using,  to  measure  ioniza¬ 
tion  rates  of  less  than  about  100/cm.  To  measure  such  low  ionization 
rates,  it  would  be  necessary  to  construct  junctions  with  breakdown 
voltage  of  the  order  of  500  -  1000  v  and  of  the  order  of  50  microns 
thick.  The  ionization  behavior  of  electrons  and  holes  would  be  of 
interest  to  predict  the  multiplication  behavior  of  high-voltage  p-n 
junctions  and  to  compare  with  theoretical  calculations. 

During  the  course  of  our  measurements,  a  minor  attempt  was  made 

to  determine  the  effect  of  crystal  imperfections  on  ionization  rate. 

Within  our  experimental  error,  the  ionization  rates  were  not  changed 

16  2 

by  neutron  bombardment  at  the  level  of  10  /cm  .  These  experiments  were 
not  sufficiently  conclusive  to  be  definitive. 

B.  MEASUREMENT  OF  ELECTRONS  EMITTED  INTO  VACUUM 

The  number  and  energy  distribution  of  electrons  emitted  into  vacuum 
from  a  reverse-biased  p-n  junction  is  an  important  clue  to  mechanisms 
of  energy  loss  by  hot  electrons.  Fairly  extensive  work  has  been  done  in 
studying  the  emission  of  hot  electrons  across  the  silicon  electron 
affinity  of  3.8  ev  into  vacuum  [Ref.  4].  This  work  has  been  analyzed 
quantitatively  to  give  a  simple  model  for  the  4  -  5-v  electron.  Such  an 
energetic  electron  has  a  lattice  collision  every  50-70  A,  giving  up 
energy  of  the  order  of  0.06  ev  (optical  phonon  energy) , and  has  a  collision 
with  the  valence  electrons  approximately  every  200  A,  resulting  in  a 
secondary  electron-hole  pair.  It  must  be  emphasized  that  the  behavior 
described  here  may  not  apply  very  well  to  the  electrons  that  cause  the 
actual  secondaries  in  a  p-n  junction.  We  have  shown  that  an  electron 
requires  a  threshold  energy  of  1.8  ±  0.1  ev  and  a  hole  requires  a 
threshold  energy  of  2.4  ±  0.1  ev  to  cause  an  ionization.  It  is  the 
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nature  of  a  threshold  phenomenon  that  Its  probability  increases  rapidly 
near  the  threshold.  Thus  the  ionization  probability  of  electrons  in  the 
range  1.7  to  2  ev  may  be  much  less  than  one  ionization  for  every  200  A 
of  travel. 

An  attempt  was  made  to  measure  the  electron  emission  and  energy 
distribution  with  a  very  much  lower  surface  barrier.  A  few  layers  of 
cesium  on  pure  silicon  can  give  a  work  function  of  the  order  of  1.5  ev. 
We  prepared  some  junctions  similar  to  the  ones  used  in  the  study  of 
vacuum  emission  across  the  silicon  work  function  and  obtained  vacuum 
emission  with  only  2  v  applied  to  the  junction.  It  was  not  possible, 
however,  to  make  the  measurements  that  are  analogous  to  the  ones  made 
with  non-cesiated  surfaces  because  the  vacuum  system  tended  to  develop 
a  shunt  leakage  path  in  the  presence  of  the  cesium  vapor.  The  lack  of 
success  does  not  mean  that  the  measurement  cannot  be  made.  The  analysis 
of  electron  emission  with  cesium  on  the  surface  is  an  experiment  that 
must  be  done  with  somewhat  more  extensive  precautions  to  avoid  leakage 
currents  than  we  exercised.  The  results  will  be  very  interesting,  as 
they  will  give  a  measure  of  ionization  probability  for  electrons  typical 
of  those  in  an  avalanching  p-n  junction. 

C.  THEORY 

Theoretical  considerations  have  been  given  for  interpreting  multi¬ 
plication  experiments,  ionization  thresholds,  ionization  rates  [Refs. 
3-5],  and  electron  emission  into  vacuum.  These  considerations  will 
not  be  repeated.  A  theoretical  difficulty  does  exist  in  relating  the 
theories  of  nonlinear  mobility  and  the  theories  of  ionization.  These 
theories  start  at  opposite  ends  of  the  electric-field  range  and  must 
be  extrapolated  toward  each  other.  If  the  extrapolations  were  within  an 
order  of  magnitude  in  predicted  electron  temperature,  we  could  be  quite 
content  that  the  phenomenon  of  hot  electrons  was  well  understood. 
Unfortunately,  the  work  of  Conwell  [Ref.  6],  and  of  Stratton  [Ref.  7], 
who  start  at  the  low-field  end,  predicts  electron  temperatures  of  the 
order  of  2000  to  3000  °K  at  electric  fields  of  2x10^  v/cm.  These 

electron  temperatures  are  not  obtained  in  the  ionization  theories  of 

5  , 

Shockley  or  Wolff  until  one  reaches  electric  field  of  about  2x10  v/cm. 
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The  electron  temperature  varies  approximately  as  the  square  of  electric 
field  in  the  high-field  range. 

A  brief  description  of  the  physical  processes  at  the  two  ends  of 
electric-field  range  will  be  given.  In  the  Conwell  and  Stratton  approach, 
it  is  supposed  that  electrons  interact  with  long-wavelength  acoustic 
phonons  and  with  optical  phonons.  The  optical -phonon  energy  is  about 
0.063  ev  in  silicon.  The  mobility  associated  with  long-wavelength  , 
acoustic-mode  scattering  has  a  T  3^2  temperature  dependence.  The 
mobility  associated  with  optical  scattering  has  a  much  steeper  dependence, 
since  the  optical  phonon  energy  exceeds  the  value  of  kT/q  in  the  range 
of  interest.  The  actual  temperature  dependence  is  fitted  by  taking  a 
mixture  of  acoustic  and  optical  phonon  scattering.  This  mixture  is  then 
used,  as  illustrated  in  the  appendix,  to  calculate  the  field  dependence 
of  velocity.  As  the  carriers  heat  up,  the  probability  of  opt ic a 1 -phonon 
interaction  increases.  A  good  fit  to  the  velocity  vs  electric  field  is 
obtained  up  to  fields  of  several  thousand  volts  per  centimeter. 

In  the  ionization  theory,  it  is  assumed  that  all  of  the  energy  loss 
to  the  lattice  is  through  optical -phonon  collisions,  and  the  mean  free 
path  for  optical -phonon  emission  is  adjusted  to  give  the  right  amount  of 
ionization.  A  value  of  70  A  must  be  taken  to  fit  the  ionization  data. 

This  value  results  in  much  too  low  an  electron  temperature  to  be  in 
agreement  with  the  extrapolated  low-field  values  even  though  it  is  in 
reasonable  agreement  with  the  spectrum  of  light  emitted  from  avalanching 
p-n  junctions. 

There  are  several  possible  ways  of  resolving  this  discrepancy.  As 
the  electrons  become  more  energetic,  new  scattering  processes  can  set  in. 
Thus,  in  silicon  the  <111>  minima  are  about  0.39  ev  above  the  <10Q> 
minima  so  that  the  density  of  end  states  is  much  greater  for  the  more 
energetic  electrons.  It  is  also  possible  that  the  assumption  of 
Maxwellian  distribution  is  not  valid  at  the  very  high  fields. 

Finally,  it  is  necessary  to  consider  the  possibility  of  scattering 
by  short -wavelength,  transverse,  acoustical  phonons.  There  is  a  high 
density  of  states  of  such  phonons,  and  their  possible  contribution  to 
electron  scattering  has  been  ignored.  This  last  possibility  is  the  one 
that  is  most  easily  tested.  The  appendix  gives  details  of  a  calculation, 
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following  the  method  of  Stratton,  which  includes  the  transverse  acoustic 
phonon.  The  result  of  this  inclusion  is  that  it  is  possible  to  obtain 
a  better  fit  to  the  velocity-vs-f ield  curve  than  if  such  scattering  is 
ignored.  The  electron  temperature  is  reduced  somewhat  but  is  still  too 
high.  Our  conclusion  must  be  that  further  theoretical  work  is  necessary 
to  resolve  this  problem. 
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APPENDIX  A.  THEORETICAL  VARIATION  OF  ELECTRON  MOBILITY  WITH 
EFFECTIVE  ELECTRON  TEMPERATURE  AND  APPLIED  ELECTRIC  FIELD  IN 

SEMICONDUCTORS 


1.  Introduction 

The  theoretical  variation  of  electron  mobility  with  the  applied 
field  in  semiconductors  has  been  studied  by  several  authors  [Refs.  7,  8]. 
Shockley  has  calculated  the  variation  of  mobility  with  the  applied  field 
in  covalent  semiconductors  assuming  the  lattice  scattering  is  mainly  due 
to  the  long-wavelength,  acoustic  phonons.  His  result  is  in  reasonable 
agreement  with  room-temperature  experiments  on  n-type  germanium  and 
silicon  only  up  to  103  v/cm.  This  fact  is  verified  by  the  experimental 
results  given  by  Ryder  [Ref.  9],  Gunn  [Ref.  10],  and  Gibson  [Ref.  ll] . 
Stratton  has  applied  Frflhlich's  and  Paranjape's  theory  [Refs.  12,  13], 
assuming  that  the  inter-electronic  collisions  essentially  determine  the 
distribution  function,  to  derive  the  variation  of  electron  mobility  with 
applied  field.  In  his  formulations,  he  completely  neglected  the  inter¬ 
valley  scattering  mechanism  due  to  energetic  phonons  near  the  zone  edge. 
Application  of  his  result  to  n-type  germanium  is  in  good  agreement  with 
experimental  data  up  to  3.5x10s  v/cm.  It  is  believed,  however,  that  the 
corresponding  electron  temperature  T  is  much  higher  than  the  actual 
situation. 

The  purpose  of  this  report  is  to  include  the  effect  of  the  inter¬ 
valley  scattering  mechanism  due  to  energetic  phonons,  both  acoustic  and 
optical.  Ohm's  law  for  electronic  conductors  is  strictly  obeyed  only  in 
the  limit  of  vanishing  applied  electric  field,  when  the  electron  gas  and 
lattice  vibrations  are  in  quasi-thermal  equilibrium.  Application  of  an 
electric  field  raises  the  mean  kinetic  energy  of  the  electrons,  and  the 
rate  at  which  energy  is  supplied  to  the  electrons  by  the  field  must  be 
balanced  by  the  rate  at  which  energy  is  transferred  from  the  electrons 
to  the  lattice  by  electron-phonon  collisions.  This  is  always  true 
before  avalanche  breakdown  occurs.  If  the  increase  of  the  mean  kinetic 
energy  is  sufficiently  large  to  make  the  electron  gas  become  warm  or 
even  hot,  the  electron  mobility  will  vary  with  the  applied  field  and 
appreciable  deviation  from  Ohm's  law  is  expected  to  occur.  The  study 
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of  a  field-dependent  electron  temperature  has  been  carried  out  by 

Frdhlich  and  Paranjape  in  considerable  detail  [Refs.  12,  13].  In  the 

latter  paper,  Frdhlich  and  Paranjape  have  pointed  out  that,  if  the 

density  of  electrons  is  sufficiently  high,  one  can  assume  that  the 

energy  and  momentum  transfers  are  controlled  by  inter-electronic 

collisions.  The  distribution  function  of  the  electrons  in  momentum 

space  will  then  be  in  internal  equilibrium  at  an  effective  electron 

temperature  T  that  depends  on  the  applied  field  F  and  is  greater 

than  the  lattice  temperature  T  .  The  whole  distribution  drifts  in  the 

direction  of  the  applied  field  with  a  velocity  equal  to  p  /m, 

O' 

We  assume  that  the  inter-electronic  collisions  essentially  determine 
the  electronic  distribution  function  f(p).  Then  the  distribution  must 
be  a  Maxwellian  distribution  that  may  be  displaced  in  momentum  space, 
i.e.  , 

I  P  P  |  2 

f(p)  =  A  exp^-  . j,  Jf(p)d3p  =  N,  (1.1A) 

where  the  normalization  constant  A  is  independent  of  p  but  is 
determined  by  the  total  number  of  electrons  N.  With  the  form  of  the 
distribution  function  assumed,  A  is  evaluated  to  be  N/(2jtmkT)3/2 . 

The  distribution  function  f(p)  contains  two  parameters  pQ  and  T, 
which  will  be  determined  from  the  rate  of  momentum  and  energy  transfer 
to  the  lattice.  If  f(p)  is  developed  in  spherical  harmonics  with  the 
axis  in  the  direction  of  the  applied  field,  then  it  can  be  expanded  as 

f(p)  =  fQ(p)  +  fj(p)  Pj^cos  e)  +  ...  ( 1 . 2A ) 


where  0  is  the  angle  between  p  and  the  applied  field  F,  and 


fl(p)  =  "P0  -£ 


(1.3A) 


(1.4A) 
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when 


2mkT 


«  1. 


(1.5A) 


Equation  (l.5A)  implies  that  the  drift  velocity  is  much  less  than 
the  random  thermal  velocity.  The  current  density  "j  and  the  drift 
velocity  v  are  related  to  pq  by 

"j  =  Nev  =  Ne(po/m),  (l.6A) 


which,  in  turn,  can  be  derived  as  functions  of  F. 

2 .  Calculation  Assuming  Inter-Electronic  Collisions  Predominate 

The  distribution  function  f(p)  in  momentum  space  must  satisfy  the 
kinetic  equation 


(2.1A) 


where  the  subscripts  F,  L  and  e  respectively  denote  the  rates  of 
change  of  f(p)  due  to  an  external  field,  to  collisions  with  the 
lattice  vibrations  including  both  intra-valley  and  inter-valley  scat¬ 
tering  mechanisms,  and  to  the  inter-electronic  collisions.  In  steady- 
state 


=  0 ♦  (  2 , 2A ) 

Since  inter-electronic  collisions  conserve  energy  and  momentum,  it 
follows  that 

^  5  (if),  -  “o  £  EP  (if),  -  0  (2-3A) 

P  P 
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where 


£ 

P 


From  Eqs .  (2.1A)  through  (2.3A),  we  obtain 


I 


(2.4A) 


and 


(2 . 5A) 


Assuming  the  distribution  function  f(p)  is  continuous  in  momentum 
space,  it  has  been  shown  [Ref.  13]  that 


(2.6A) 


and 


eFp  N 
o 

m 


(2.7A) 


where  the  field  F  is  in  the  z  direction.  Equation  (2.6A)  is  related 
to  steady-state  momentum  of  the  electrons  and  Eq.  (2.7A)  is  the  rate  of 
supplying  energy  from  the  electric  field.  Frdhlich  and  Paranjape  have 
also  considered  the  detailed  processes  of  absorption  and  emission  of 
phonons  by  electrons  to  calculate  (3f/5t)^.  The  result  is 

(!)l3  X[g°(p)i+P°Pl  (C°S  8)  gl(p)iL  (2'8A) 


SEL-64-009 


-  22 


where 


hiph  =  - 


go(p)i 


V  fo(p) 


-  —  J  f  (p)  f  Q  dq  B(q)  N 
2ntf4  P  °  J  11 

l-exp(o>o-tD) ]  +  +  [exp  cuo~exp  o>]_}  (2.9A) 


2nti* 


C  Cr  /  2  mk9.v 


exp  (cd^-od) 


exp  coq-  (1- 


>■(' 


mk0. 


2p2 


exp  go 


(2.10A) 


and 

o>  =  0./T  ,  =  0,/T 

o  i'  o  i' 

where 


Tq  =  lattice  temperature 
T  =  electron  temperature. 


The  plus  sign  (+)  outside  the  bracket  indicates  the  term  con¬ 
tributed  from  absorption  of  a  phonon,  and  the  minus  sign  (-)  emission 
of  a  phonon.  The  subscript  i  corresponds  to  each  mode  of  the  vibra¬ 
tional  phonons:  B(q).  is  the  square  of  the  interaction  matrix  element 
t  h  ^ 

due  to  i  mode,  and  the  number  of  phonons  corresponding  to  different 
modes  is 


Ni  "  i^pfe  )--i ■  ®  <«>  -*“(q>  *  k6i-  (2 . ha) 

o' 

Here  q/h  is  the  wave  number  of  the  phonons  that  scatter  the  electrons. 
From  the  vibrational  spectra  of  covalent  semiconductors  Si  and  Ge 
[Refs.  14,  15]  we  expect  that  each  of  these  modes  will  contribute  both  in 
the  intra-valley  and  inter-valley  scattering  mechanisms.  From  Eqs.  (2.4A) 
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through  (2.8A)  and  making  use  of  Eq.  (l.6A) 


we  obtain 


JF  = 


gQ(p)id3p 


(2.12A) 


.  intra 


df 


inter 


St/ta  Wio^aJ 


bt\ 


inter" 


d3p 


i 


r  r^f\intra  /bi\intra  [bt\inter 

J  pz  \St/a  +TOo  +  \3t/ta  + 


(bt\inter 

\5t/io,ia_ 


d3p 


(2.13A) 


Here  we  use  the  property  of  orthogonality  relation  for  Legendre 
polynomials;  only  the  zero-order  term  Po(cos  0)  contributes  to  the 
energy  integral,  while  the  first-order  term  P^cos  0)  contributes 
only  to  the  momentum  integral.  The  superscripts  stand  for  different 
modes  of  the  phonons. 

The  upper  and  lower  limits  of  integration  for  q  in  Eqs .  (2.9A) 
and  (2.10A)  are  found  to  be  [Ref.  13] 

^±11  =  P  +  (p2±2mk9^)^  u(p2±2mk9^)  (2.14A) 

q+/0  =  P  -  (p2±2mk9^)^  U(p2±2mk0^)  (2.15A) 


for  the  bracket  marked  +  and  -  respectively.  u(x)  is  equivalent  to  a 
unit  step  function  that  is  zero  when  x  is  negative  and  unity  when  x 
is  positive.  The  0^  are  the  equivalent  Debye  temperatures  correspond¬ 
ing  to  different  modes  of  the  phonon  both  at  zone  center  and  near  the 
zone  edge  in  the  reduced  Brillouin  zone  representation. 
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For  intra-valley  scattering  of  low-energy,  long-wavelength,  acoustic 
phonons,  the  upper  and  lower  integration  limits  of  Eqs.  (2.9A)  and 
(2.10A)  should  be  2p  and  0  respectively.  Since  we  know  from  the 
vibrational  spectra  and  the  order  of  magnitude  of  the  mean  electron 
energy  that  qQ  >  2p,  after  simple  calculation  we  have 


,  > intra  fo^  4p2  (  1  To  \ 

P  a  ='  p  4mkTQ4a  \,kT  ’  kX  J 


(2.16A) 


gx(p) 


intra 

a 


fQ(p) 

mkT 


32mp3 

5k2T 


(1- 

J_V  4p if 

\T 

2V  3k2 

( 2 . 17A) 


Here  we  make  use  of  the  relations  [Refs.  13,  16] 


and 


B(q)intra  =  _2C^q__ 

3  9M's  V  n. 

1 


e (q)  =  qs  ^ 


qQ  =  4i(6nzni)1//3 
4iM  :S2q® 

l  =| - 2 

a  8n  C2m2kT 


/  (2.18A) 


V 

2ith4 


mB(q) 


intra 

a 


_ . 

"  4 mkT  I 
o  a 


Now,  from  Eqs.  (2.16A)  and  (2.17A),  and  applying  the  relation  of 
(2.8A)  into  Eqs.  (2.12A)  and  (2.13A),  we  obtain 


(2.19A) 
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In  calculating  the  result  of  Eq.  (2.20A),  we  retain  only  up  to  the 
first  order  of  s  in  Eq.  (2.17A).  The  ratio  of  the  second  term  to  the 
first  term  in  the  right-hand  side  of  Eq.  (2.20A)  is  of  the  order  of 
10“4  at  room  temperature  for  covalent  semiconductors  Si  and  Ge.  Thus 
the  second  term  is  completely  negligible. 

In  order  to  calculate  the  effects  that  come  from  inter-valley 
scattering  near  the  zone  edges,  we  assume  that  the  scattering  mechanisms 
due  to  the  energetic  phonons  either  at  the  zone  center  or  on  the  zone 
edge  are  of  the  same  form.  Then  we  may  substitute  the  upper  and  lower 
limits  of  integration  from  Eqs .  (2.14A)  and  (2.15A)  into  Eqs .  (2.9A)  and 
(2.10A),  also  making  use  of  the  relations  [Refs.  16,  17] 


to  obtain 
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go(p)i 


f  (p)  2E 
O  p 

p  (2Ni+l) 


P 

mkT 


f0(p) 

\(e) 


He) 

T^Je) 


(2.22A) 

(2.23A) 


where 


1 


(2E/m)'/s  ^[^(kS./E)]^  (n.+1)  [  l-(k9.  /e)  ]'*  U[l-(k8./E)] 
l.  2N.  +  1 

1  l 


(2.24A) 


Again  carrying  out  the  energy  and  the  momentum  integrals  using 
Eqs.  (2.12A)  and  (2.13A),  we  obtain 


"I EP  (sf)i d3p  =  Ne  h  ^r.'('2N.+i)  (f)s  [(Ni+1)  e_a)-Ni]  JG(“) 


(2.25A) 


®i d3p  =  Nev  M  (£)  {[(vi}  e_“+Ni]  Ji(cw) 


+  a>  JQ(co)  (N^+l)  e  j  (2.26A) 


where 

4e^i 

=  - r 7  • 

3(2jimkT  r 
'  o' 

J1(cu)  and  Jq(co)  are  two  integrals  that  can  be  evaluated  in  terms 
of  modified  Bessel  functions  of  the  second  kind*  as 


We  would  like  to  thank  Professor  G.  E.  Latta  who  helped  us  evaluate 
the  integrals  in  exact  forms. 
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Jo(co)  =  J  e“x[x(x-Ku)]/a  dx  =  |  exp(|)  K^) 

(2.27A) 


00 

J^cn)  =  ^  e"x[x3(x-Ko)]^  dx 
o 


=  a,exp(|){Kl(|)-  f  [Kl(|)-Ko(f)]) 


Now,  combine  all  the  individual  terms  in  Eqs .  (2.12A)  and  (2.13A) 
together,  lumping  the  terms  due  to  optical  phonons  into  one  term  with 
a  weighting  factor  of  1.7.  Since  0^  =  0.80 ,  detailed  analysis  shows 
that  this  factor  is  quite  reasonable.  Actually,  the  weighting  factor 
here  is  not  important  at  all  since  there  are  still  some  adjustable 
parameters  in  the  final  results  that  should  be  determined  from  the 
temperature  dependence  of  the  measured  electron-mobility  data.  However, 
it  does  indicate  the  physical  picture  implied  in  the  formulations.  The 
final  results  are 


vF  = 


.  i_  k9a  / jrU 

3*  ^;(2Na+D  \To) 


t(Na+l) 


e"0-  N  ]  J  (a) 
a  o  ' 


o  1.7  k0 
8 _ o 

+  3n  mp.  (2N  +l) 


(f 

'  n 1 


[(No+l) 


"7 


Nq]  J0(7) 


(2.28A) 


32  /  T  16 


v  =  9npa  (to/  +  9n  l2Na 


TIT^  iff  ([(Na+1>  V0)  +  aJo^) 

,ra  '  n' 


(V1*  *  i  isrib;  (ff  (“V') *  "o1 


+  yJQ(y)  (N0+l)  e“7)  .  (2.29A) 
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The  field  dependence  of  the  drift  velocity  v  can  be  evaluated 
from  these  results  if  we  can  eliminate  the  effective  electron  temperature 
T  between  Eqs.  (2.28A)  and  (2.29A).  Three  special  cases  of  zero,  very 
weak,  and  very  strong  applied  fields,  will  be  considered  in  more  detail: 

a .  For  the  Zero-Field  Case 

In  the  limit  of  vanishing  applied  field,  the  mobility  is  given  by 


— =  lim  (-) 
^cA  F->0  V 


(2.30A) 


Substituting  the  relations  of  Eq .  (2.31A)  back  into  Eqs.  (2.28A)  and 
(2.29A)  leads  to 


3k0  .  .J4 

(2.32A) 
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f  -  IT-  (f  f  +  ^  (f f  ([(Na+l)e^Na]Jl(a)+aJo(a)(Vi)e^} 

^aA  'o'  ^aA  '  o' 

/4 

['.IWWH'1  *  5^  (f) 

♦7Jo(7)(No«)e")'){2No[Ji(>'0)»(ro/2)^0(ro)]r1  .  (2.33A) 

The  variation  of  the  zero  field  mobility  p  with  the  lattice 
temperature  will  be  considered  as  follows.  Since  &  is  inversely 
proportional  to  T  ,  pa^  can  therefore  be  written  as 

UaA  =  b(7o)3/2  (2.34A) 


where 


b 


3n  e  1 

8  03/2  (2nmk)y2 


constant . 


Furthermore,  from  Eqs.  (2.21A),  also  using  the  relations 
✓HK/2  a  qQ  and  qQS  3  k@ ,  where  0  is  the  corresponding  Debye  tempera¬ 
ture  of  the  low-energy  acoustic  phonons  near  zone  center,  we  may  derive 
the  following  relations 


K  i  ^yDe2 

i*  =  4r  To9a 
a 

K  x  (2No+l)02 


( 2 . 35A) 


(2.36A) 


where  r  =  C2/l8D2  is  a  dimensionless  constant  of  the  order  of  a  few 
hundredths  [Refs.  16,  18].  Therefore,  we  have,  from  Eq .  (2.31A), 
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(2.37A) 


,/  4r6?  i 

-  '•W  — 

aS 


,,  4r6‘ 

--  '  TT 


{N0[ji(7o)+(7o/2)J0(7o-)]}“1  (2.38A) 


where 


If  we  introduce  the  reduced  quantities 


cA 
b  ’ 


R1  " 


4r  02 


R  = 

2  4r92 


(  2 . 39A) 


and  substitute  them  into  Eq.(2.30A),  we  obtain 


i  =  (y  )-a/2  + 

M  K7o' 


aR2[  Jx( a70) +( a70/2) JQ(  a70) ] 

(70)/s(ea7o-l) 

1 • 7Rj[ Jj( 7q)+( 70/2) Jq( 7q) 1 

(y)He7o-l) 


(2.40A) 


The  plots  of  M  vs  1/7  for  two  extreme  cases  are  given  in  Figs. 
10  a  and  b.  One  extreme  case,  which  is  =  0  and  R^  adjustable,  is 
equivalent  to  Herring's  result  [Ref.  19], 

b.  For  the  Very  Weak  Field  Case 

We  may  assume  that  T  =  Tq  +  AT,  and  expanding  Eqs.  (2.32A)  and 
(2.33A)  in  power  series  of  S  =  AT/Tq  «  1,  retaining  up  to  the  first 
order  of  8,  we  thus  eliminate  8  between  the  two  equations.  We  then 
solve  v  in  terms  of  F.  The  result  gives 


H  = 


(2.41A) 
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iG.  10.  PLOT  OF  THE  REDUCED  ZERO-FIELD  MOBILITY- TEMPERATURE 
CURVES  ACCORDING  TO  EQ.  (2.40A). 


where 


.  f  3k9  ii  .  a  J  (a  ) 

F?  =  -i-  As=n— i  - 2-^-2- 

1  2  \  \  ( /r»\ 


VW  )  r  FjA  r  “oJ1fao>«0fao-‘)J0(°o)-1 

2»  ^oA  ^  L  *  V'Vt(<V2>Jol“o'> J 

’ilr0M?j2Kb')  -]]  ■  <2'42A) 


(2.42A) 


c .  For  a  Very  Strong  Field  Case 

If  the  applied  field  is  sufficiently  strong  to  make  T  »  9q  and 

T  »  0  ,  i.e.,  y  «  1,  a  «  l.then  Eqs.  (2.28A)  and  (2.29A)  reduce  to 
a 


■f[Sr(# 


,/a  1.7kd 


3mp^(2Na+l)  \To'  3mti 


•  KBo  (_T.fl 
o(2No+1>  \To /  _ 


( 2 . 43A) 


(i  •  *  •  a  (if 


(2.44A) 


The  best  fit  of  the  calculation  to  the  experimental  data  of  n-type 
silicon,  suggests  that  the  third  term  in  the  right-hand  side  of  Eq. 
(2.43A)  dominates  almost  all  the  way  through  the  applied  field  range 
before  breakdown  occurs.  (The  same  also  occurs  for  n-type  germanium, 
which  is  not  presented  in  this  report).  Therefore,  we  are  unable  to  get 
F^  dependence  of  the  drift  velocity  v  at  high  field,  which  is  true 
only  when  the  first  term  in  the  right-hand  side  of  Eq.  (2.43A)  pre¬ 
dominates.  We  do  get  the  result  that  the  theoretical  fitted  curve 
almost  comes  to  saturation  at  high  field  range. 

3 .  Application  of  the  Results  to  n-Type  Silicon 

The  band  structure  of  silicon  based  on  a  quantum-theoretical 
calculation,  has  been  successfully  carried  out  by  F.  Herman  [Ref.  20] 
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and  D,  P.  Jenkins  [Ref.  21].  The  conduction  band  of  Si  has  six 
equivalent  minima,  which  occur  at  points  lying  in  the  <100>  directions 

—P  —P  **♦> 

and  having  values  of  K  equal  to  about  0.8  K  ,  where  K  is  the 

max  max 

*-> 

value  of  K  corresponding  to  the  zone  boundary  in  these  directions. 

The  effective  mass  of  electrons  in  the  conduction  band,  which  is  actually 
associated  with  variation,  is  assumed  constant  in  the  fitting  process. 

Its  value  is  calculated  from  the  relation 


1  =  1  /2_ 
m  3  \mt 


(3.1A) 


Substituting  the  values  nr  =  0.19  m  ,  m»  =  0.98  m  [Ref.  22],  we 

t  o  l  o 

obtain  m  =  0.26  mQ.  The  sound  velocity  in  n-type  silicon  can  be  evaluated 
by  comparing  the  vibrational  spectra  of  Ge  and  Si  [Refs. 14, 15]  and  also 
by  making  use  of  the  sound  velocity  of  5.4x10  cm/sec  in  Ge  [Refs.  8,  23] 
to  give  s  =  7.05x10s  cm/sec.  The  effective  Debye  temperatures  corre¬ 
sponding  to  the  optical  phonons  at  zone  center  and  transverse  acoustic 
phonons  near  zone  edge  can  be  evaluated  from  the  vibrational  spectra 

[Ref.  14],  The  results  are  0  =  690  °K  and  0  =  200  °K.  The  experi- 

o  a 

mental  result  of  n-type  silicon  given  by  Ludwig  [Ref.  24]  will  be  used 
to  fit  the  temperature  dependence  of  the  zero  field  mobility.  Applying 
their  result  yields 

PcA  =  2.1x10®  T"2-5  cm2/v-sec.  (3.2A) 


If  150  °K  5  T  S  350  °K,  graphical  analysis  shows  that  T-25 
°  o 

dependence  can  be  fitted  to  within  ±10  percent  with  the  best  fit  of  the 
field  dependence  of  drift  velocity  v.  The  final  result  is  shown  in 
Fig.  11  with  R^  =  1.5  and  R^  =  1.5  that  leads  to 


3/2 


H aA  *  3.16xl03  cm2 /v- sec 

(690  V*  exp(200/To)  -  1 

\  T0 )  J1(200/Tq)  +  (I00/To)  Jo(200/Tq)  Cm  /V_sec 


( 3 . 3A) 


MaA  =  7-26xl°3 


|i  .  =  1.24x10s 
oA 


’moV‘  exp(690/TQ)  -  1 

,  V  •Y690/V  *  <-34S/T0)  J0(6®°/Tc,>  “  /v"sec' 
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T.  <°K> 


FIG.  11.  FITTED,  TEMPERATURE-DEPENDENT  CURVE 
OF  ZERO-FIELD  ELECTRON  MOBILITY  TO  n-TYPE 
SILICON  WITH  b  =  3.16xl03  cm2/v-sec. 
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Using  the  values  of  p^  ,  p^,  and  deduced  in  Eq.  (3.3A), 

it  is  possible  to  evaluate  F2  from  Eq.  (2.42A).  If  we  assume  that 

the  fitting  of  the  temperature-dependent  curve  at  zero  field  can  be 

extrapolated  down  to  100  °K,  then  we  may  calculate  the  quadratic 

departure  of  electron  mobility  due  to  the  applied  field  for  both 

T  =  99  °K  and  T  =  300  °K.  The  results  are  listed  in  Table  1,  and 
o  o 

compared  to  Brown’s  experimental  result [Ref.  26], 


TABLE  1.  EVALUATION  OF  F2  AT  VARIOUS  LATTICE 

1 

TEMPERATURES 


T 

^aA 

^aA 

Ft 

(v/cm)2 

(  °k) 

^’aA 

^oA 

Theoretical 

Experimental 

99 

2.06 

0.216 

1.51x10s 

— 

300 

2.91 

3.18 

3.92xl07 

1.85xl07 

We  can  also  evaluate  the  theoretical  variation  of  drift  velocity 
v  with  the  applied  field  F  by  inserting  the  deduced  values  of 
p  ,  p  ’  and  p  .  into  Eqs.  (2.32A)  and  (2.33A).  The  best  result 

aA  flA  On 

is  shown  in  Fig.  12  together  with  the  experimental  curve  obtained  by 
A,  C.  Prior  [Ref.  25],  The  result  shows  that  the  calculated,  effective, 
electron  temperature  T  is  about  six  times  higher  than  the  estimated 
actual  values  at  high  field.  The  discrepancy  between  the  theoretical 
evaluated  values  of  v  and  the  measured  data  is  about  1.5. 


TABLE  2.  THE  ESTIMATED  MEAN  ELECTRON  ENERGY  AND  EFFECTIVE 
ELECTRON  TEMPERATURE  ASSOCIATED  WITH  THE  POINTS 
MARKED  ON  THE  THEORETICAL  CURVE  IN  FIG.  12 


Number  of 
Point 

Mean  Electron  Energy 
Estimated  (ev) 

Effective  Electron 
Temperature  (°K) 

1 

0.0267 

310 

2 

0.0298 

345 

3 

0.0414 

480 

4 

0.0595 

690 

5 

0.147 

1700 

6 

0.298 

3450 

7 

0.595 

6900 

8 

1.490 

17250 
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ELECTRIC  FIELD  E  (v/cm) 


FIG.  12.  THE  VARIATION  OF  DRIFT  VELOCITY  WITH  THE  APPLIED  FIELD  FOR 
n-TYPE  Si  ACCORDING  TO  THEORY,  TOGETHER  WITH  THE  EXPERIMENTAL 
RESULT  GIVEN  BY  PRIOR.  Solid  curve  is  Prior's  data  for  4.8  ft-cm 
sample;  broken  curve  is  the  best  fit  for  the  theoretical  result. 
The  estimated  mean  electron  energy  associated  with  the  points 
marked  is  given  in  Table  2. 
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4.  Discussion  and  Conclusion 


The  result  which  shows  both  discrepancies  in  the  theoretical  calcu¬ 
lated  electron  temperature  and  drift  velocity  is  not  surprising.  Since 
we  oversimplify  the  physical  picture  implied  in  our  formulations,  there 
are  at  least  some  points  which  would  influence  the  result  to  some  extent. 
They  are: 

1.  The  form  of  distribution  function:  We  assume  the  distribution 
function  as  the  form  of  Eq,  (l.lA)  that  is  a  Maxwellian  distribu¬ 
tion  displaced  in  momentum  space.  We  expand  it  in  terms  of 
spherical  harmonics  as  shown  in  Eq.  (l.2A)  and  cut  off  right  after 
the  second  term.  Obviously,  it  is  oversimplified,  especially  for 
the  case  when  the  applied  field  is  sufficiently  high,  while  the 
distribution  function  may  be  far  distorted  from  the  quasi¬ 
equilibrium  position.  It  is  believed  that  this  is  extremely 
serious  in  the  direction  of  the  applied  field,  thus  making  the 
approximations  used  in  the  calculation  breakdown. 

2.  The  effective  mass  of  the  electrons  in  the  conduction  band:  We 
assume  the  effective  mass  of  the  electrons  in  the  conduction  band 
to  be  constant  through  the  calculation.  Actually,  this  is  not  the 
case.  The  effective  mass  of  the  electrons  in  the  conduction  band 
is  given  by  the  relation 


-  =  —  V2  E  .  (4.1A) 

m  A2  k 


According  to  the  calculated  band  structure  of  silicon,  it  is 
clearly  seen  that  the  effective  mass  of  the  electron  is  actually 
a  variable  in  the  conduction  band.  This  will  certainly  have  some 
influence  in  the  calculation. 

3.  The  scattering  mechanism  of  the  energetic  phonons  near  zone  edge: 
The  scattering  mechanism  of  energetic  phonons  near  zone  edge  is 
assumed  to  be  of  the  same  form  as  that  of  the  energetic  phonons 
at  zone  center.  It  seems  quite  reasonable  to  assume  in  this  way 
as  to  distinguish  from  the  different  scattering  mechanism  due  to 
the  low  energy,  long-wavelength  acoustic  phonons.  This  assumption 
also  makes  it  possible  to  evaluate  BCq)^,  square  of  the  inter¬ 
action  matrix.  However,  the  theoretical  justification  is  opened 
to  question  up  to  the  present.  If  more  information  can  be  obtained 
from  this  kind  of  scattering  mechanism,  we  may  be  able  to  establish 
a  clearer  physical  model  implied  in  the  formulations. 
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4.  The  scattering  mechanism  involved  between  different  minima  in  the 
conduction  band:  This  is  completely  neglected  in  our  formulation. 
However,  if  the  applied  field  is  sufficiently  high  as  to  make  the 
electrons  in  the  conduction  band  become  hot,  this  kind  of  scat¬ 
tering  may  be  able  to  occur.  The  energy  minimum  in  <111>  direc¬ 
tion  is  about  0.39  ev  higher  than  that  in  <100>  direction  in  the 
conduction  band  of  silicon,  which  is  probably  too  high  for  this 
kind  of  scattering  to  occur  even  at  high  fields.  (However,  it  is 
quite  likely  to  occur  in  germanium  since  the  energy  minima  in  the 
<000  and  <100  directions  are  only  about  0.1  ev  and  0.18  ev  higher 
than  that  in  the  <111>  direction  [Ref.  20].) 

5.  Complication  of  band  structure:  The  complication  of  band  structure 
may  result  in  different  effective  electron  temperatures  at  the 
minima  of  different  valleys  and  thus  influence  the  final  result. 

The  fitting  of  the  theoretical  calculation  to  the  experimental 
result  and  the  evaluation  of  the  quadratic  departure  of  electron  mobil¬ 
ity  with  the  applied  field  shows  that  the  scattering  mechanism  of  the 
optical  phonons  is  of  great  importance  in  the  energy  transfer  processes 
from  the  electrons  to  the  lattice,  even  when  the  mean  electron  energy 
is  less  than  that  of  the  optical  phonons.  This  may  be  explained  by  the 
reason  that  even  though  the  number  of  excited  optical  phonons  is  small, 
the  energy  transfer  due  to  these  excited  optical  phonons  is  much  greater 
than  that  due  to  the  low  energy  acoustic  phonons. 

Application  of  the  theoretical  calculation  to  n-type  silicon 
shows  that  the  results  are  still  a  factor  of  6  higher  in  electron 
temperature  and  a  factor  of  1.5  higher  in  drift  velocity  than  the 
accepted  experimental  results.  (The  result  is  even  worse  in  n-type 
germanium,  corresponding  factors  are  10  and  2  respectively. )  This 
suggests  that  the  distribution  function  and  the  scattering  mechanism  have 
to  be  considered  in  more  detail. 

However,  comparison  of  the  result  obtained  here  to  Stratton's 
(n-type  Ge)  clearly  indicates  that  the  method  applied  in  our  formulation 
results  in  almost  saturation  at  the  high  field  range;  a  result  which  is 
impossible  to  obtain  if  Stratton's  method  is  applied.  (The  improvement 
in  effective  electron  temperature  is  less  than  a  factor  of  2  in 
n-t  ype  Ge . ) 
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